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In this paper we present an algorithm for computing Grobner bases of linear ideals in a difference polynomial 
ring over a ground difference field. The input difference polynomials generating the ideal are also assumed to 
be linear. The algorithm is an adaptation to difference ideals of our polynomial algorithm based on Janet-like 
reductions. 



1. INTRODUCTION 

Being invented 40 years ago by Buchberger pQ 
for algorithmic solving of the membership prob- 
lem in the theory of polynomial ideals, the 
Grobner bases method has become a power- 
ful universal algorithmic tool for solving various 
mathematical problems arising in science and en- 
gineering. 

Though overwhelming majority of the Grobner 
bases applications is still found in commutative 
polynomial algebra, over the last decade a sub- 
stantial progress has also been achieved in appli- 
cation of Grobner bases to noncommutative poly- 
nomial algebra, to algebra of differential oper- 
ators and to linear partial differential equations 
(see, for example, book 2). As to the difference 
algebra, i.e. algebra of difference polynomials, in 
spite of its conceptual algorithmic similarity to 
differential algebra, only a few efforts have been 
done to extend the theory of Grobner bases to 
difference algebra and to exploit their algorith- 
mic power |3l4j . 

Recently, two promising applications of differ- 
ence Grobner bases were revealed: generation 
of difference schemes for numerical solving of 
PDEs |5I6| and reduction of multiloop Fcynman 
integrals to the minimal set of basis integrals [J] . 

In this note we describe an algorithm (Sect. 4) 
for constructing Grobner bases for linear differ- 
ence systems that is an adaptation of our polyno- 
mial algorithm 8 to linear difference ideals. We 
construct a Grobner basis in its Janet-like form 



(Sect. 3), since this approach has shown its com- 
putational efficiency in the polynomial case [EE]- 
We briefly outline these efficiency issues in Sect. 5. 
The difference form of the algorithm exploits 
some basic notions and concepts of difference al- 
gebra (Sect. 2) as well as the definition of Janet- 
like Grobner bases and Janet-like reductions to- 
gether with the algorithmic characterization of 
Janet-like bases (Sect. 3). We conclude in Sect. 6. 

2. ELEMENTS OF DIFFERENCE AL- 
GEBRA 

Let {y 1 , . . . , y m } be the set of indetermi- 
nates, for example, functions of n— variables 
{xi, . . . , x n }, and Q\, . . . , 9 n be the set of mutu- 
ally commuting difference operators (differences), 
e.g., 9i o yi =y3(xi,...,Xi + l,.. .,x n ). 

A difference ring R with differences 9\ , . . . , 6 n 
is a commutative ring R such that V/, g G R, 1 < 
i, j < n 

9i9j = OjOu e i o(f + g) = e i of + 9 i o g , 

i °(fg) = (9 i of)(9 i og). 

Similarly, one defines a difference field. 

Let IK be a difference field, and ffi. := 
Kjy 1 , . . . , y m } be the difference ring of polyno- 
mials over K in variables 
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Hereafter, we denote by the set of linear poly- 
nomials in R and use the notations: 

= {0" | ^eZ| }, de gi (^oy fc )= Mi , 

deg(^oj, fe ) = lMl = Er=i/^- 
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A difference ideal is an ideal I C K closed 
under the action of any operator from 0. If 
F := {fx, . . . , f k } C K is a finite set, then the 
smallest difference ideal containing F will be de- 
noted by Id(F). If for an ideal I there is F C K L 
such that / = Id(-F), then / is a linear difference 
ideal. 

A total ordering >- on the set of ^ o j/ 3 is a 
ranking if k 7 /x, v the following hold: 

9 t 9» oyi y O^oyi , 

0» o yl y 9 V o y k 9 t 9^ o yi y 9 t 9 v o y k . 

If /i y v =^> 9^ o y 3 y 9 V o y k the ranking is 
orderly. If j y k => 9^ o y 3 y 9 U o y k the ranking 
is elimination. 

Given a ranking y, a linear polynomial / G 
Ml \ {0} has the leading term a$ o y 3 , ■§ G 0, 
where $ o y 3 is maximal w.r.t. >- among all 9^ o 
y k which appear with nonzero coefficient in /. 
lc(/) := a G K \ {0} is the leading coefficient and 
lm(/) := o y 3 is the leading monomial. 

A ranking acts in Ri as a monomial order. 
If F C R L \ {0}, lm(F) will denote the set of 
the leading monomials and lnx,- (F ) will denote its 
subset for indeterminate y 3 . Thus, 

lm(F) = U7 =1 lm J (F). 

3. JANET-LIKE GROBNER BASES 

Given a nonzero linear difference ideal / = 
Id(G) and a ranking y, the ideal generating set 
G = {g%, ■ ■ ■ , g s } C Ml is a Grobner basis |2I4| of 
/ if V/ G Jflli \ {0}: 

3geG,6ee : lm(/) = 9 o lm(g) . (1) 

It follows that f E I \ {0} is reducible modulo G: 

f->f':=f-\c(f)9o(g/\ C (g)), f G I. 
a 

If /' ^ 0, then it is again reducible modulo G, 
and, by repeating the reduction, in finitely many 
steps we obtain 




can be reduced to an irreducible polynomial h, 
which is said to be in the normal form modulo F 
(denotation: h — NF(h, F)). 

In our algorithmic construction of Grobner 
bases we shall use a restricted set of reductions 
called Janet-like (cf. [Hj) and defined as follows. 

For a finite set F C Ml and a ranking y, we 
partition every lmk(F) into groups labeled by 
d ,...,di G Z> , (0 < i < n). Here [0]fc := 
lmk(F) and for i > the group [do,...,di\k is 
defined as 

{u G lm fe (F) | d = 0, dj = deg^u), 1 < j < i}. 
Denote by hi{u, hxik(F)) the nonnegative integer 

maxldeg^u) | u,v G [do, d;_i] fc } - deg;(u). 

If hi(u,lm.k(F)) > 0, then 0f< such that 

s,; := min{deg i (w) - deg^u) | 

u, v G [d , di-i]fc, deg^u) > deg^u)} 

is called a difference power for f £ F with 
lm(/) = u. 

Let DP(f, F) be the set of difference powers 
for / G F, and F) := 9 \ 6 be the subset of 
with 

9:={0"\30 u eDP(f,F) : M -^gZ| }. 

A Grobner basis G of / = Id(G) is called Janet- 
like HJ if V/ G inRi \{0}: 

3. 9 gG,0g J(g,G) : lm(/) = 9 o lm(g) . (2) 

This implies J— reductions and the J — normal 
form NFj(f,F). It is clear that condition (J2J) 
implies Note, however, that the converse is 
generally not true. Therefore, not every Grobner 
basis is Janet-like. 

The properties of a Janet-like basis are very 
similar to those of a Janet basis 9 , but the for- 
mer is generally more compact than the latter. 
More preciously, let GB be a reduced Grobner 
basis [2], JB be a minimal Janet basis, and JLB 
be a minimal Janet-like basis of the same ideal for 
the same ranking. Then their cardinalities satisfy 



Similarly, a nonzero polynomial h G Ml, whose 
terms are reducible (if any) modulo a set F C Ml, 



Card(GB) < C&rd(JLB) < C&rd{JB), (3) 
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where Card abbreviates cardinality, that is, the 
number of elements. 

Whereas the algorithmic characterization of a 
Grobner basis is zero redundancy of all its S— 
polynomials |1I2| . the algorithmic characteriza- 
tion of a Janet-like basis G is the following con- 
dition (cf. (HI): 

Vg G G, d e DP(g, G) : NFj(#og, G) = . (4) 

This condition is at the root of the algorithmic 
construction of Janet-like bases as described in 
the next section. 

4. ALGORITHM 



This algorithm is an adaptation of the polyno- 
mial version [H] to linear difference ideals. It out- 
puts a minimal Janet-like Grobner basis which 
(if monic, that is, normalized by division of each 



polynomial by its leading coefficient) is uniquely 
defined by the input set F and ranking >~. Cor- 
rectness and termination of the algorithm follow 
from the proof given in [8]; in so doing the dis- 
placement of some elements of the intermediate 
sets G into Q at step 13 provides minimality of 
the output basis. The algorithm terminates when 
the set Q becomes empty in accordance with 

The subalgorithm Normal Form(p, G, y) per- 
forms the Janet-like reductions (Sect. 3) of the in- 
put difference polynomial p modulo the set G and 
outputs the Janet-like normal form of p. As long 
as the intermediate difference polynomial h has a 
term Janet-like reducible modulo G, the elemen- 
tary reduction of this term is done at step 4. As 
usually in the Grobner bases techniques [2], the 
reduction is terminated in finitely many steps due 
to the properties of the ranking (Sect. 2). 



Algorithm: Normal Form(p, G, y) 

Input: p£li\ {0}, a polynomial; GcMi\ 

{0}, a finite set; y, a ranking 
Output: h = NFj[p, G), the ./-normal 
form of p modulo G 
1: h := p 

2: while h ^ and h has a monomial u with 
coefficient b G K J— reducible modulo G 
do 

3: take j € G such that u = 6 1 o Im(<?) 

with 6~t G J(lm( 5 ),lm(G)) 
4: h:=h/b-0fo( g /lt(g)) 
5: od 

6: return h 



An improved version of the above algorithm can 
easily be derived from the one for the involutive 
algorithm (jj] if one replaces the input involutive 
division by a Janet-like monomial division |SJ and 
then translates the algorithm into linear differ- 
ence algebra. In particular, the improved version 
includes Buchberger's criteria adjusted to Janet- 
like division and avoids the repeated prolonga- 
tions 8 13 o g at step 16 of the algorithm. 



Algorithm: Janet-like Grobner Basis^, y) 



Input: F C R^\{0}, a finite set; y, a ranking 


Output: G, a Janet-like basis of Id(_F) 


1 


choose / G F with the lowest lm(/) 




w.r.t. y 


2 


G-= {/} 


3 


Q:=F\G 


4 


do 


5 


h := 


6 


while and /i = 0do 


7 


choose p G Q with the lowest lm(p) 




w.r.t. y 


8 


Q:=Q\ M 


9 


h := Normal Form(p, G, >-) 


10 


od 


11 


if ft, ^ then 


12 


for all g G G such that lm(g) = 9^ o 




lmf», ImI > do 


13 


Q:=QU{g}: G:=G\{g} 


14 


od 


15 


G := GU{h} 


16 


Q := Q U { 9? o g | g G G, 9? G 




DP(g,G)} 


17 


fi 


18 


od while Q ^ 


19 


return G 
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5. COMPUTATIONAL ASPECTS 

The polynomial version of algorithm Janet- 
like Grobner Basis implemented in its im- 
proved form in C++ (Hj has disclosed its high 
computational efficiency for the standard set of 
benchmarks 2 . If one compares this algorithm 
with the involutive one [H] specialized in Janet 
division, then all the computational merits of the 
latter algorithm are retained, namely: 

• Automatic avoidance of some useless reduc- 
tions. 

• Weakened role of the criteria: even with- 
out applying any criteria the algorithm is 
reasonably fast. By contrast, Buchberger's 
algorithm without applying the criteria be- 
comes unpractical even for rather small 
problems. 

• Smooth growth of intermediate coefficients. 

• Fast search of a polynomial reductor which 
provides an elementary Janet-like reduction 
of the given term. It should be noted that 
as well as in the involutive algorithm such 
a reductor, if it exists, is unique. The fast 
search is based on the special data struc- 
tures called Janet trees [§]. 

• Natural and effective parallelism. 

Though one needs intensive benchmarking for lin- 
ear difference systems, we have solid grounds to 
believe that the above listed computational mer- 
its hold also for the difference case. 

As this takes place, computation of a Janet- 
like basis is more efficient than computation of a 
Janet basis by the involutive algorithm f§j- The 
inequality 10) for monic bases is a consequence of 
the inclusion [Hj: 

GB C J LB C JB . (5) 

There are many systems for which the cardinality 
of a Janet-like basis is much closer to that of the 
reduced Grobner basis than the cardinality of a 
Janet basis. Certain binomial ideals called toric 



form an important class of such problems. Toric 
ideals arise in a number of problems of algebraic 
geometry and closely related to integer program- 
ming. For this class of ideals the cardinality of 
Janet bases is typically much larger than that of 
reduced Grobner bases [B]. For illustrative pur- 
poses consider a difference analogue of the simple 
toric ideal |8lllj generated in the ring of differ- 
ence operators by the following set: 




The reduced Grobner basis for the degree-reverse- 
lexicographic ranking with 9 X >~ 9 y y 9 Z >- 9 W is 
given by 

{ ^\ — ^ y ^z,8%9 w — 9 y ,9 x 9y — 9 z 9 Wl 9y — 9 x 9 z 9 2 w } . 

The Janet-like basis computed by the above al- 
gorithm contains one more element 9 X 9 W — 9 y 
whereas the Janet basis adds another six extra 
elements to the Janet- like basis [B]. 

The presence of extra elements in a Janet ba- 
sis in comparison with a Janet-like basis is ob- 
tained because of certain additional algebraic op- 
erations. That is why the computation of a Janet- 
like basis is more efficient than the computation 
of a Janet basis. Both bases, however, contain the 
reduced Grobner basis as the internally fixed ||] 
subset of the output basis 3 . Hence, having any of 
the bases computed, the reduced Grobner basis is 
easily extracted without any extra computational 
costs. 

6. CONCLUSION 

The above presented algorithm is implemented, 
in its improved form, as a Maple package 
and already applied to generation of difference 
schemes for PDEs and to reduction of some loop 
Feynman integrals (see some examples in |10|). 
The last problem for more than 3 internal lines 
with masses is computationally hard for the cur- 
rent version of the package. 

One reason for this is that the Maple imple- 
mentation does not support Janet trees since 
Maple does not provide efficient data structures 
for trees. 



2 See Web page http : //in vo . j inr ■ ru| 



3 In the improved versions of the algorithms. 
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Another reason is that in the improved version 
of the algorithm there is still some freedom in 
the selection strategy for elements in Q to be 
reduced modulo G. Though our algorithms are 
much less sensitive to the selection strategy than 
Buchberger's algorithm, the running time still 
depends substantially on the selection strategy: 
mainly because of dependence of the intermedi- 
ate coefficients growth on the selection strategy. 
To find a heuristically good selection strategy one 
needs to do intensive benchmarking with differ- 
ence systems. In turn, this requires an exten- 
sive data base of various benchmarks that, unlike 
polynomial benchmarks, up to now is missing for 
difference systems. 

In addition to our further research on improve- 
ments of the Maple package, we are going to im- 
plement the difference algorithm in C++ as a 
module of the open source software GINV avail- 
able on the Web site http : //invo . j inr . ru 

The comparison of implementations of poly- 
nomial involutive algorithms for Janet bases in 
Maple and in C++ ^3 shows that the C++ code 
is of two or three order faster than its Maple coun- 
terpart. Together with efficient parallelization of 
the algorithm this gives a real hope for its prac- 
tical applicability to problems of current interest 
in reduction of loop integrals. 
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